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Abstract. In this article, we study the axisymmetric surface diffusion flow 
(ASD), a fourth-order geometric evolution law. In particular, we prove that 
ASD generates a real analytic semiflow in the space of (2 + «)-little-H61der reg- 
ular surfaces of revolution embedded in R 3 and satisfying periodic boundary 
conditions. We also give conditions for global existence of solutions and prove 
that solutions are real analytic in time and space. Further, we investigate the 
geometric properties of solutions to ASD. Utilizing a connection to axisym- 
metric surfaces with constant mean curvature, we characterize the equilibria 
of ASD. Then, focusing on the family of cylinders, we establish results regard- 
ing stability, instability and bifurcation behavior, with the radius acting as a 
bifurcation parameter for the problem. 



1. Introduction 

The central focus of this article is the development of an analytic setting and 
rigorous results for the axisymmetric surface diffusion flow (ASD) with periodic 
boundary conditions. We establish well-posedness of ASD and investigate geometric 
properties of solutions, including characterizing equilibria and establishing their 
stability, instability and bifurcation behavior. We establish and take full advantage 
of maximal regularity for ASD. Most notably, with maximal regularity we gain 
access to the implicit function theorem, a very powerful tool in nonlinear analysis 
and dynamical systems theory. We begin with a motivation and derivation of the 
general surface diffusion flow, of which ASD is a special case, and we introduce the 
general outline of the paper. 

The mathematical equations modeling surface diffusion go back to a paper by 
Mullins [35] from the 1950s, who was in turn motivated by earlier work of Herring 
[24] . Both of these authors investigate phenomena witnessed in sintering processes, 
a method by which objects are created by heating powdered material to a high 
temperature, while remaining below the boiling point of the particular substance. 
When the applied temperature reaches a critical point, the atoms on the surfaces of 
individual particles will diffuse across to other particles, fusing the powder together 
into one solid object. In response to gradients of the chemical potential along the 
surface of this newly formed object, the surface atoms may undergo diffusive mass 
transport on the surface of the object, attempting to reduce the surface free energy. 
Given the right conditions - temperature, pressure, grain size, sample size, etc. - 
the mass flux due to this chemical potential will dominate the dynamics on the 
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surface, and it is the resulting morphological evolution of the surface which the 
surface diffusion flow aims to model. We also note that the surface diffusion flow 
has been used to model the motion of surfaces in other physical processes (e.g. 
growth of crystals and nano-structures) . The article [H] contains the formulation of 
the model which we present below, which is set in a more general framework than 
the original model developed by Mullins. 

f.f . The Surface Diffusion Flow. From a mathematical perspective, the gov- 
erning equation for motion via surface diffusion can be expressed for hypersurfaces 
in arbitrary space dimensions. In particular, let T C K™ be a closed, compact, 
immersed, oriented Riemmanian manifold with codimension 1. Then we denote by 
H = H(T) the (normalized) mean curvature on F, which is simply the sum of the 
principle curvatures on the hypersurface, and Ar denotes the Laplace-Beltrami op- 
erator, or surface Laplacian, on T. The motion of the surface T by surface diffusion 
is then governed by the equation 

V = A r H, 

where V denotes the normal velocity of the surface T. A solution to the surface 
diffusion problem on the interval J C K+, with € J, is a family {T(t) : t 6 J} of 
closed, compact, immersed hypersurfaces in K™ which satisfy the equation 



(1.1) 



(V(r(t)) = A m H(T(t)), teJ:=J\ {0}, 

\r(o) = r 0) 



for a given initial hypersurface Tq. It can be shown that solutions to (1.1) are 



volume-preserving, in the sense that the signed volume of the region Q enclosed 



by the surface T is preserved along solutions. Additionally, (1.1) is surface-area- 
reducing. 

Well-posedness, for any dimension n > 2, was established by Escher, Mayer and 
Simonett in |20) . where it was also shown that the (n — l)-dimcnsional spheres 
are asymptotically stable equilibria. Meanwhile, Mayer and Simonett |34j demon- 
strate the existence of initially embedded hypersurfaces which are driven to self- 
intersection under the surface diffusion flow. Beyond these initial results, the lit- 



erature lacks general analytic results regarding the behavior of solutions to (1.1), 
e.g. occurrence of singularities and conditions under which solutions breakdown in 
finite time. An important feature of the surface diffusion flow that Escher, Mayer 
and Simonett exploit in order to obtain well-posedness results is the fact that the 
equation has a quasilinear structure and the linear part of the equation exhibits 
maximal regularity properties on appropriately chosen function spaces. These fea- 
tures will also play an important role in our analysis of ASD, for which we establish 
a robust theory for local and global well-posedness. 



1.2. Axisymmetric Surface Diffusion (ASD). For the remainder of the paper, 
we will focus our attention on a special case of the surface diffusion flow. Namely, 
we consider the case of T C R 3 an embedded surface which is symmetric about 
an axis of rotation (which we take to be the a;-axis, without loss of generality) 
and satisfies prescribed periodic boundary conditions on some fixed interval L of 
periodicity (we take L = [— 7r,7r] and enforce 2tt periodicity, without significant loss 
of generality). In particular, the axisymmetric surface T is characterized by the 
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parametrization 

T = i(x,r(x)cos(6),r(x)sm(6)) : xeR, 6 e [— tt, tt]|, 

where the function r : R — > (0, oo) is the profile function for the surface T. 
Conversely, a profile function r : M — > (0, oo) generates an axisymmetric surface 
r = r(r) via the parametrization given above. 

Utilizing the explicit parametrization for axisymmetric surfaces, we can recast 
the surface diffusion problem as an evolution equation for the profile functions 
r = r(t). In particular, one can see that the surface T(r) inherits the Riemannian 
metric 

g = (1 + r 2 x ) dxAdx + r 2 dO A d9, 
from the embedding r c — > R 3 , with respect to the surface coordinates (x, 6); where 
the subscript f Xi := d Xi f indicates the derivative of / with respect to the indicated 
variable Xi. It follows that the (normalized) mean curvature of the surface is T-L{r) = 
K\ + «2, where 



Hi 



and H2 — 



(l + r|)3/2 



are the azimuthal and axial principle curvatures, respectively, on T(r). Meanwhile, 
the Laplace-Beltrami operator on F and the normal velocity of T = T(t) are 



Ar( r ) 



I 



ryfl + r 
rt 



d x 



X 



do 



v(t) = 

Finally, plugging these terms into the equation ( |1.1[ ) and simplifying, we arrive at 
the expression 



(1.2) 



1 



d x 



d x 



r(t,x + 2tt) = r(t,x), t > 0, x e M, 

^ r(0, x) = r (x), iel, 

for the periodic axisymmetric surface diffusion problem. To simplify notation in 
the sequel, we define the operator 



t > 0, x e 



(1.3) 



1 



G(r) := - a 



d x H(r) 



which is formally equivalent to the right hand side of the first equation in (1.2 1 



The first investigations of evolution of an axisymmetric surface via surface diffu- 
sion can be traced back to the work of Nichols and Mullins [36l|37] in 1965, where 
one can already see some of the benefits of this special case of the surface diffu- 
sion problem. In particular, Mullins and Nichols are able to take advantage of the 
symmetry of the problem in order to develop an adequate scheme for numerical 
techniques. Moreover, Mullins and Nichols are already predicting the finite-time 
pinch-off of tube like surfaces via surface diffusion flow, a feature similar to the 
mean curvature flow and a natural phenomenon to study in exactly this axisym- 
metric setting. Following this seminal work by Nichols and Mullins, there are many 
publications investigating ASD. Researchers continued to study pinch-off behavior 



4 



JEREMY LECRONE AND GIERI SIMONETT 



using numerical methods, c.f. [TOl [TTJ [12l |T5l [311 E2] , developing schemes for the 
continuation of solutions after the change of topology that occurs at the moment of 
pinch-off. Meanwhile, many researchers have focused on the numerical investiga- 
tion of stability/instability and bifurcation behavior of cylinders, which are natural 
equilibria of ASD, under perturbations of various types, c.f. [TOl [12] , though 
we note that these investigations lack a rigorous analytic framework. For more 
background on ASD see the paper [7] by Bernoff, Bertozzi, and Witelski, which 
has served as a motivation for the overall nature and scope of results which have 
been justified for ASD, via numerical and analytic techniques, which we can hope 
to establish in our functional-analytic setting. 

In Section[2j we prove existence of solutions to ( 1.2 1, which, to the best of the au- 
thors' knowledge, is the first analytic well-posedness result in the literature for the 
axisymmetric surface diffusion flow with periodic boundary conditions. In particu- 
lar, we establish existence and uniqueness of maximal solutions which are analytic 
in time and space for positive time, with a prescribed singularity at time t = 0, for 
initial conditions which are (2 + a)-little-H61der continuous in space. Additionally, 
we establish conditions for global existence and regularity of the semiflow induced 
by (1.2 1. With these well-posedness results established, all of which depend heav- 
ily upon the theory developed in [29 and the wcll-poscdncss results for quasilincar 
equations with maximal regularity proved by Clement and Simonett [5], we then 
go on to investigate more general dynamic properties of solutions. In Section [3] we 
characterize the equilibria of ASD using a result of Delaunay [17] and Kenmotsu 
[26] . which describes all of the constant mean curvature surfaces in the axisymmet- 
ric setting. In Section [4j we prove that the family of cylinders with radius r* > 1 
are asymptotically, exponentially stable under a large class of nonlinear perturba- 
tions, which maintain the same axis of symmetry and satisfy the prescribed periodic 
boundary conditions. To prove this stability result, we linearize the equation for 
ASD and see that the spectrum of the linearized operator is contained in the left 
half of the complex plane, although the spectrum will always contain zero as an 
eigenvalue. Hence, we reduce the equation, by essentially eliminating non-volumc- 
preserving perturbations, in order to eliminate the zero eigenvalue. Then we use 
results regarding maximal regularity on exponentially weighted function spaces to 
generate the desired exponential stability for the reduced equation, which is then 
transferred back to the (full) ASD problem via a lifting operator. 

We proceed in the remaining sections with results regarding the instability of 
cylinders with radius < r+ < 1 and the existence of branches of bifurcating 
equilibria which intersect the family of cylinders at radii r+ = \/t, for every I € 
N := {1,2,...}. The instability result of Section [5] makes use of a contradiction 
technique reminiscent of results from the theory of ordinary differential equations, 
c.f. Priiss and Wilke [13] • More precisely, by isolating the linearization of the 
governing equation, one takes advantage of a spectral gap and associated spectral 
projections in order to derive necessary conditions for stable initial data, which in 
turn lead to a contradiction. We also refer to Priiss, Simonett and Zacher [42] and 
Priiss, Simonett and Wilke [40] for related results. Finally, in Section [6] we apply 
classic methods of Crandall and Rabinowitz [13 to establish bifurcation results 
for ASD. However, even in the reduced setting mentioned above, developed while 
proving stability results, we find that the eigenvalues associated with the problem 
are not simple. We restrict our attention to surfaces which are even (symmetric 
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about the surface x — 0) and satisfy prescribed regularity and periodicity, similar 
to a method used by Escher and Matioc [IS]. In this setting, the problem has 
simple eigenvalues and we are able to derive bifurcation results and apply them 
back to the full ASD problem via a posteriori symmetries of equilibria. We note 
that, rather than using the restriction to even functions in order to apply the 
results of Crandall and Rabinowitz |13j . we could also have chosen to apply more 
general bifurcation techniques, such as the methods contained in the manuscript of 
Kiclhofer [57J Section 1.19], in order to generate (higher dimensional) bifurcation 
results. 

1.3. Maximal Regularity. One tool that we will use extensively throughout the 
paper is the property of (continuous) maximal regularity, also called optimal reg- 
ularity in the literature. Maximal regularity has received a lot of attention in 
connection with parabolic partial differential equations and evolution laws, c.f. 
[21 El SI El EH [3U1 EH SH [S] . Although maximal regularity can be developed in a 
more general setting, we will focus on the setting of continuous maximal regularity 
and direct the interested reader to the references [H [30] for a general development 
of the theory. 

Let n G (0, 1], J := [0,T], for some T > 0, and let E be a (real or complex) 
Banach space. Following the notation of [9], we define spaces of continuous functions 
on J := J \ {0} with prescribed singularity at 0. Namely, define 

BUd-pi^E) := jit G C{j,E) : [t ^ t 1 '^)] G BUC(J,E) and 
( L4 ) Urn f 1 "'>(i)|| B = o} 1 (0,1) 

IM| Bl _ ti i-sup^MOiu, 

teJ 

where BUC denotes the space consisting of bounded, uniformly continuous func- 
tions. It is easy to verify that BUCx-^J, E) is a Banach space when equipped with 
the norm || • ||bi_ • Moreover, we define the subspace 

BUCl_^(J,E) := {uGC^j.E) : u, u G BUCx-^ J, £)} , A 1 <= (0,1) 
and we set 

BUC {J,E) := BUC{J,E) BUCq{J,E) := BUC l {J,E). 

Now, if Ei and Eq are a pair of Banach spaces such that E\ is continuously 
embedded in Eq, denoted E% ^ Eq, we set 

E Q (J):= BUd-^Eo), A*e (0,1], 

Ei (J) := BUCl^J, E ) n BUd.^J, E,), 

where Ei(J) is a Banach space with the norm 

IM| El(J) ^suptWHOii^ + M^ii^Y 

te J 

It follows that the trace operator 7 : Ei(J) — > Eq, defined by -/v :— v(0), is well- 
defined and we denote by 7E1 the image of 7 in E , which is itself a Banach space 
when equipped with the norm 

||a;|| T Ki := inf I IMKi;./) : v G Ei (J) and = x\. 
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For a bounded linear operator B G C{E\,Eq) which is closed as an operator 
on E , we say (Eo(J),Ei(J)J is a pair of maximal regularity for B and write B G 
if 

+ B ' 7 ) G Asom(El(J),Eo(J) X 7 El), 

where Ci Som denotes the space of bounded linear isomorphisms. In particular, 
(E (J),Ei(J)) is a pair of maximal regularity for B if and only if for every (/, u Q ) G 
Eo( J) x 7E1, there exists a unique solution u G Ei( J) to the inhomogeneous Cauchy 
problem 

u(t) + Bu(t) = f(t), tej, 
u(0) = u . 

Moreover, in the current setting, it follows that 7E1 = (Eo, ^ i.e. the trace 
space 7E1 is topologically equivalent to the noted continuous interpolation spaces 
of Da Prato and Grisvard, c.f. EE El ESS- 



2. Well-Posedness of (1.2 1 



When considering the surface diffusion problem, the underlying Banach spaces 
Eq and E\ in the formulation of maximal regularity will be spacial regularity classes 
which describe the properties of the profile functions r(t). We proceed by defining 
these regularity classes. We define the one-dimensional torus T := [— tt,w], where 
the points — ir and w are identified, which is equipped with the topology generated 
by the metric 

dr(x, y) := va\n{\x - y\, 2ir - \x - y\}, x, y £ T. 

There is a natural equivalence between functions defined on T and 27r-periodic func- 
tions on ]R which preserves properties of (Holder) continuity and differentiability. 
In particular, we will be working with the so-called periodic little-Holder spaces 
/i CT (T), for g G R + \ Z. Definitions and basic properties of periodic little-Holder 
spaces, as well as details on the connection between spaces of functions on T and 
27r-periodic functions on K can be found in and the references therein. For the 
readers convenience, we provide a brief definition of /i CT (T) below. 

For k e N := NU{0}, denote by C fc (T) the Banach space of /c-times continuously 
diffcrcntiable functions / : T — > R, equipped with the norm 

k k 

:( 



\f\\c H r) ■■= E H/ (j) llc(T) := £ ( sup \f^(x)\ 
i=o j=o 



Moreover, for a G (0,1) and k G No, we define the space C k+a (T) to be those 
functions / G C k (T) such that the a- Holder seminorm 



[/«] , := SUP V""<*J-f'yi 
x, v <ET a T (x,y) 



is finite. It follows that C k+a (T) is a Banach space when equipped with the norm 

\\f\\o>+-m :=ll/llc fc( T) + [/ (fe) ] 
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Finally, we define the periodic little-Holder space 

I, " ! ' ' ! T ):={f€ C k+a (T) : lim sup l/ W (^) ~ / W fa)l _ q v 

0<d r (x,y)<8 

for k G No and a G (0, 1) which is a Banach algebra with pointwise multiplication 
of functions and equipped with the norm |j • H^M-a := || • ||c fe + a (T) inherited from 
C k+a (T). For equivalent definitions and more properties of the periodic little- 
Holder spaces, see [351 Section 1]. 



In order to make explicit the quasilinear structure of (1.2 1, we reformulate the 
problem. By expanding the governing equation we arrive at the formally equivalent 
problem 



(2.1) 



\d t r(t,x) + [A(r(t))r(t)](x) = f(r(t,x)), t > 0, x G T, 
\r(0,x) =r (x), x G T, 

where, for appropriately chosen functions p, 



(2-2) A( P ) := — * dt + Wl + ^Z^) # 

(I + Pl) 2 p(l + P lf 



is a fourth-order differential operator with variable coefficients over T and 
/ on ,/ s Pi - 1 - 1 2 3 - \bpj 3 



p2(l + /5 2 )2 — p(l +/9 2)3^ (1 + ^)4 p3 (1 + p 2) 

is a R- valued function over T. Looking at these formal expressions, one can deduce 
several properties that the functions p must satisfy in order to get good mapping 
properties for / and A. In particular, we want to choose p such that p(x) ^ for 
all x € T, also we want that the spacial derivatives p x and p xx make sense and the 
products p 2 , p 3 , pp x , etc. have desired regularity properties. With these conditions 
in mind, we proceed with our well-posedness result. 

2.1. Existence and Uniqueness of Solutions. Fix a G (0,1) and define the 
spaces of R- valued little-Holder continuous functions 

(2.4) E :=h a (T), Ei := /i 4+Q (T), and E, := (E , E^, 

where (•, -)° ^ for p G (0, 1), denotes the continuous interpolation functor of Da 
Prato and Grisvard, c.f. [14] or [2]. It is well-known that the little-Holder spaces 
are stable under this interpolation method, in particular we know that 

E^ = h^ L+a {T) (up to equivalent norms), for Ap + a £ Z, 

c.f. [291 EQ]. Further, let V be the set of functions r : T ->• R such that r(x) > 
for all x G T and define := V (1 E^ for p G [0, 1]. We note that is an open 
subset of Ep, for all p G [0, 1]. 

Lemma 2.1. Let p G [1/2,1]. Then 

{A,f)eC u (v^, MR v (E 1 ,E )xE \ for v e (0,1], 
where C u denotes the space of real analytic mappings between Banach spaces. 
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Proof. Fix fi £ [1/2, 1] as indicated. This result relies on the fact that the little- 
Holder spaces /i CT (T) are Banach algebras with pointwise multiplication of functions. 

CLAIM 1: A(p) E MR V (E 1 ,E ) for p G V M , v E (0, 1]. This claim will follow 
from , though the setting of that paper differs slightly from the current setting 
and warrants some discussion. First, for p G define the coefficients 

1 , , , x 2p z (l + pi - ippxx) 



M/?) '■= 77~ — Tv? and & s(p) 



t^ + Pl? p(l + pl) 3 ' 

so that A(p) — 64 (p) + b 3 (p)d^.. By our choice of p, it follows that C 
h 2+a (T, E), so that 64, 63 E Eo and A(p) is a uniformly elliptic differential operator. 
By [53 Theorem 5.2] we conclude that 

A(p) eMR u (h 4+a (T,C),h a (T,C)), ve (0,1], 

where we utilize the notation h k+a (T, C) to be clear that the space consists of de- 
valued functions over T, and does not coincide with the spaces E^ being considered 
herein. However, h k+a (T, C) does coincide with the complexification of h k+a (T, M.) 
(up to equivalent norms) and it is a straightforward exercise to see that the property 
of maximal regularity continues to hold under restriction to the subspaces /i CT (T, R). 

CLAIM 2: The operation of inversion is real analytic from Vq into Eq, i.e. 

T z : [r4l/r] E C"(V ,E ). 

Fix tq E Vq and choose a > so that tq(x) > a for all ieT. If r E Eo is chosen 
so that || r — ro\\c(T) < a i then the representation 

111 / ro(x) — r(x) s 



r(x) r (x)(l+ r W;$ x) ) r o( x )fo\ 

holds for x eT arbitrary, where the last equality follows by an elementary geometric 
series argument. Hence, the given power series represents the function 1/r point- 
wise for x eT. Moreover, it follows from the algebraic structure of Eq that 

JL ( r _ r )n 

^ : =E n+l ^0, fOT fceN - 

n=0 r 



Finally, if ||ro — t\\e < a/2, then 
(r - r) n 



71=0 uiiEq n =0 



||^||n+l - 2v „n+l ~ Z / a > 



< > ^ ^° < 

n=0 

which demonstrates the power series converges absolutely in the topology of Eq and 
the claim follows. 



CLAIM 3: The operations of differentiation and multiplication are real analytic 
in the setting of little-Holder spaces, i.e. 

d x : [r^r x }EC"(h° +1 {T),h°(T)), a E R+ \ Z, 

T m : \{r, s) ^ rs] E C"(E x Eo, Eq). 

This claim follows since d x and T m are bounded, linear and bilinear (respectively) 
on the indicated spaces. 

The regularity of the mappings A and / now follows from the preceding claims, 
with the additional observation that the mapping A : V M — > £(Ei,E ) inherits the 
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regularity of the coefficients 03 , 04 : — > Eq and the fact that MR^Ei, E$) is an 
open subset of C(E , Ex), c.f. [HI Lemma 2.5(a)]. □ 



With Lemma |2.1| established, we can take full advantage of the well-posedness 
results for quasilinear parabolic equations presented in the article [9] of Clement 
and Simonett. In particular, we conclude results regarding local existence and 
uniqueness of solutions with continuous dependence on initial data, as well as max- 
imal solutions and conditions for global existence. However, we give only a limited 
presentation, focusing on those results which will be of most direct use to us in the 
sequel, and we refer the interested reader to [5] for further details on well-posedness 
of ( |2~T] ). 

Before we can properly state a result on maximal solutions, we need to introduce 
one more space of functions from an interval J C R+ to a Banach space E, with 
prescribed singularity at zero. Namely, if J = [0, a) for a > 0, i.e. J is a right-open 
interval containing 0, then we set 

Ci-^E) := {u G C(j,E) : u G BUCi-^{[0, T],E), T < sup J}, 

Cl.^E) := {u G C\j,E) Ci_ M (J,£)}, M G (0, 1], 

which we equip with the natural Frechet topologies induced by BUCi-^dO, T], E) 
and BUC^^T^E), respectively. 

Theorem 2.2 (Existence and Uniqueness). Fix a G (0, 1) and take fi G [1/2, 1] so 
that 4/i + a ^ 7L. For each initial value ro G := h i ^ +a (T) n [r > 0], there exists 
a unique maximal solution 

r(;r ) G Cl_ M (J(r ), h a (T)) n Ci_ M (J(r ), /i 4+Q (T)), 

where J(ro) = [0,^ + (ro)) M + denotes the maximal interval of existence for initial 
data ro . Further, it follows that 

V~ (J J(r Q )x{r } 

r £V M 

is open m R + x and 95 : [(t,ro) > r(f,r"o)] is an analytic semiflow on V^, i.e. 
using the notation tp (ro) := v?(i, ro), ifte mapping cp satisfies the conditions 

• tp G C(T>, Vfi) 

• w° = idv 

(2 5) 

• <fi s+t (r Q ) = ip l o tp s (r ) for < s < t+(r ) and0<t< t+(ip s (r Q )) 

• (p(t, •) G C u (V t , 7 P ) /or i G R+ with V t := {r G V„ : (t, r) G V} ^ 0. 



Proof. In case /i G [1/2,1), the result follows from Lemma 2.1 and [3J Theorems 
4.1, 5.1 and 6.1]. When /.i = 1 we note that the existence and uniqueness of a 
maximal solution 

r(-, r ) G C\ J(r ), £ ) n C(J(r ), Si) 
follows from [9J Theorem 4.1(b)]. However, for the semiflow properties, we will 
consider (1.2) as a fully nonlinear equation, and apply results of Angenent 5J. In 
particular, for r G V\ we use the representation G(r) = — A(r)r + f(r) and (2.2)- 
(2.3) to see that the Frechet derivative DG has the structure 

DG(r) = -TTT"^ dt + ]T B k (r) d k x , 
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where the coefficients Bk{r) 6 Eq for every r 6 V\, k = 0, . . . , 3. From this 
computation it foliows that —DG(r) is a uniformly elliptic operator from E\ to 
Eq and so, using the results of 2!)j as in Claim 1 of Lemma 2.1 above, we see 
that -DG{r) G MRi{Ex,Eq) for all r G Vi. Now the fact that (|1.2|) generates an 
analytic semiflow on Vi follows from [5] Corollary 2.9]. □ 

The results contained in [S] also give the following conditions for global solutions. 
We have separated this result from the previous existence result because breakdown 
of solutions to (2.1 1, in particular an analytic investigation of pinch-off behavior of 
certain solutions, is an open and interesting topic. 

Theorem 2.3 (Global Solutions). Let r £ := h^ +a (T) n [r > 0] for /j G 

(1/2, 1], such that 4/i + a Z, and suppose there exists < M < oo so that, for all 
t G J(r ) 

• r(t,r )(x) > l/M, Va; G T, and 

• IK*> r o)lk 4 *+«(T) < M, 

then it must hold that i + (fo) = oo, so that r(-,r ) is a global solution. Conversely, 
if ro G Vfi and t + (ro) < oo 7 i.e. the solution breaks down in finite-time, then one, 
or both, of the conditions stated must fail to hold. 

We can also state the following result regarding analyticity of the maximal so- 
lutions r(-,r ) in both space and time. 

Theorem 2.4 (Regularity of Solutions). Under the same assumptions as in The- 
orem \2.S\ it follows that 

(2.6) r(,r„)eC u ((0,(+(r„))xT) for all r G //G [1/2,1]. 

Proof. Here we rely on an idea that goes back to Masuda [53] and Angenent [S] |5] 
to introduce parameters and use the implicit function theorem to obtain regularity 
results for solutions, see also [2T1 [22l [23] . 

First, for a e R let T a : T — > T be the translation operator, where T a (x) denotes 
the unique element in T that is in the coset [x + a] G M/27rZ of (x + a). T a 
naturally acts on functions u G C(T,R) by virtue of (T a u)(x) :— u(T a (x)). As in 
[2"2"] one shows that, for a G M, the family of translations {T ta ■ t G M} induces a 
strongly continuous group of contractions on any of the spaces E^, with infinitesimal 
generator A a given by 

D(A a ) = h 1+i » +a (T,R), A a = ad x . 

Let ro G Vfj, be fixed, and let 

r = r(; r ) G Cl_Jj(r ),E ) n Ci_^(J(r ), Bi) 



be the unique solution to (2.1 1 on the maximal interval of existence J(fo) = 
[0,t+(r )). Let £i G (0,i+(r )) be fixed and set / := [0,t t ]. Then there ex- 
ists S > such that (1 + X)t G J(r"o) for all (i, A) G I X (-5,5). Finally, for 
(A, a) eW := (-S, S) 2 we set 

rx, a {t) :=T ta r({l + X)t), teJ; 

i-e. rA,a(^^) = + A)£,T ta (a:)) for (i, a;) G I X T. One verifies that 

^A,a e E x (7) := BVCl_ M (r,Eo) n BUC^i^Ej. 
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Moreover, since the nonlinear mapping [r i— > G(r)] is equivariant with respect to 
translations, i.e. TJ, G(r) = G(T{, r) for any b G R, we obtain that r\ a is a solution 
of the parameter-dependent equation 



(2.7) 



\d t v = (l + X)G(v) + ad x v, t > 0, 

\«(0) = ro, 

on the time interval /. 

Now, for U(7) := E X (J) n C(7, V) we define 

$ : U(J) x W -> E (Z) x (A, a)) = (ftu - (1 + \)G(v) - afl.w, 7« - r ) , 

where E (J) := BUCi-^I, E ), and we note that $(rA, Q! (A, a)) = (0, 0). Moreover, 



$ G C 



' (u(I) x Vt 7 , E (J) x , £>!*(r, (0, 0)) = - OG(r), 7 j 



where we use the same notation for r = r(-,ro) and its restriction to the time 
interval /. Exactly as in the proof of [H Theorem 6.1] one shows that 

Di$(r,(0,0)) G A.om(E 1 (J),E (2) x E„). 

Finally, according to the implicit function theorem, c.f. [16j Theorem 15.3] or [18j 
(10.2.1)], there exist a neighborhood of r in Ex (J) and a neighborhood of (0,0) 
in R 2 , which we will again denote by U(J) and W , respectively, and a mapping 
5 G C u (W,Ex(I)) such that 

(A, a)) = (0, 0) if and only if v = g(\, a) 

whenever (v, (A, a)) G U(J) x VF. We conclude that g(\,a) — r\^ a and 

(2.8) [(A,a)^r A ,J GC"(W,U(7)). 
For to G (0, ti) and xo G T fixed, we see that 

(2.9) [(A, a) ^ r((l + A)t , ^ a(^o))] G C"(W, M), 

and the assertion follows since (toi^o) can be chosen arbitrarily. □ 

Remark 2.5. The preceding results can be slightly weakened to allow for arbitrary 
values of /i G (1/21], i.e. without eliminating the possibility that + a G Z, 
by taking initial data from the continuous interpolation spaces (i?o,£'i)° which 
coincide with the Zygmund spaces over T. 

3. Characterizing The Equilibria of ASD 



With well-posedness of (1.2 1 established, we move on to investigate geometric 



properties of solutions. We begin our analysis of the long-time behavior of solutions 



by characterizing and describing the equilibria of ( 1.2 1. For this characterization, we 
make use of a well-known, strict Lyapunov functional for the surface diffusion flow, 
namely the surface area functional, and a characterization of surfaces of revolution 
with prescribed mean curvature, as presented by Kenmotsu [2"6l . 



Recalling the operator G, as expressed by (1.3) and taking it to be defined on 



V\ C h 4+a (T), one will see that the set of equilibria of (1.2) coincides with the 



null set of G. Although, from the well-posedness results of the previous section, 



we know that we can consider (1.2) with initial conditions in h 2+a (T), upon which 



the operator G is not defined, one immediately sees that all equilibria must be 



in h A+a (T) (in fact, by Theorem 2.4 we can even conclude that equilibria are in 
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C°°(T)). More specifically, if we define equilibria to be those elements r£ V1/2 
V n h 2+a (T), such that the maximal solution r(-,f) satisfies 



r(t,r)=r, t > 0, 

then it follows immediately that f g /i 4+Q (T) and G(f) 
characterizing the elements of the null set of G. 
Consider the functional 



0. Now, we proceed by 



S(r) := / r(x)y/l + r*(x)dx, 



which corresponds to the surface area of the generated surface T(r) and is a strict 
Lyapunov functional for (4.1 1. Indeed, if r = r(-,ro) is a solution to (1.2) on the 



interval J(fo), then (suppressing the variable of integration) 



dtS(r(t)) = 



r{t)r x (t) 



G(r(t)) dx 



r(t) 



d x H(r(t)) H(r(t)) dx 



r(t) 



T y/l + r*{t) 



(d x H(r(t))) 2 dx, 



t € J(r ) \ {0}, 



where we use integration by parts twice and eliminate boundary terms because of 
periodicity. Notice that the expression is non-positive for all times t £ J(r"o) \ {0}. 
Moreover, if f is an equilibrium of ( |1.2| it follows that d x W.(f) is identically zero on 
T. Meanwhile, notice by the definition of the operator G that G(r) — whenever 
d x H(f) — 0. Hence, we conclude that S(r) is a strict Lyapunov functional for ( 1.2 1, 



as claimed, and we also see that the equilibria of (1.2) are exactly those functions 
f £ h i+a (T) for which the mean curvature function H(f) is constant on T. 

The axisymmetric surfaces with constant mean curvature have been character- 
ized explicitly by Kenmotsu in |26j . In particular, we see that all equilibria of 



(1.2) are so-called undulary curves, and the unduloid surfaces, which are generated 



by the undulary curves by revolution about the axis of symmetry, are stationary 



solutions of the original surface diffusion problem (1.1) 



Theorem 3.1 (Delaunay [17] and Kenmotsu |26j). Any complete surface of revo- 
lution with constant mean curvature T-L is either a sphere, a catenoid, or a surface 
whose profile curve is given (up to translation along the axis of symmetry) by the 
parametric expression, parametrized by the arc-length parameter s£K, 
(3.1) 



R(s;H,B) := 



l + Bsin(Ht) 



y/l + B 2 + 2B sm(m) 



dt 



y/l + B 2 + 2Bsm(Hs) 



\H\ 



Remarks 3.2. We can immediately draw several conclusions from Theorem |3 . 1 1 and 
characterize the equilibria of (1.2). We use the notation R(H,B) to denote the 
graph in M 2 of the parametric expression R(- ; H, B). 

a) Although the curves R(H, B) are well-defined for arbitrary values B E K 
and H ^ 0, it is not difficult to see that, up to translations along the x-axis, 
we may restrict our attention to values H > and B > 0, c.f. [551 Section 
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2]. However, in the sequel we will consider the unduloids in the setting of 
even functions on T, for which we will benefit by allowing B to range over 
all of E. 

b) When |J3| = 1, R(T~L, 1) corresponds to a family of spheres controlled by the 
parameter H, The spheres are a well-known family of stable equilibria for the 
surface diffusion flow, c.f. |20j . however their profile curves are outside of our 
current setting because they fail to be continuously differentiable functions on 
all of T. Moreover, we should note that the spheres represented by R(H, ±1) 
are in fact a connected family of spheres, or a chain of pearls (see Figure lljrl 
for which even general techniques for break down, as the manifold is 
singular at the points of intersection. These families of connected spheres 
may be interesting objects to investigate in a weaker formulation of ASD, 
but they fall outside of the current setting. 

c) Catenoids, or more precisely the generating catenary curves (which are essen- 
tially just the hyperbolic cosine, up to scaling), do not fall into the current 
setting because they fail to satisfy the periodic boundary conditions, c.f. 
Figure [TJ 

d) In case \B\ > 1, the curve R(H,B) is called a nodary (see Figure [2]), which 
cannot be realized as the graph of a function over the x-axis and hence falls 
outside the current setting. 

e) For values < \B\ < 1, R(W,B) is a family of undulary curves, which 
generate the unduloid surfaces. The undulary curves are representable as 
graphs of functions over the x-axis, which are strictly positive for B in the 
given range (see Figure [3|. In fact, the case B = corresponds to the 



cylinder of radius 1/H. Hence, by Theorem 3.1 above, we conclude that all 



equilibria of (1.2) fall into the family of undulary curves. 

f) Notice that the curve R(H, B) is always periodic in both the parameter s 
and the spacial variable x. In order to ensure that the curve satisfies the 



27r-periodic boundary conditions enforced in (1.2) (which we emphasize is a 
condition regarding periodicity over the variable x and not the arc-length 
parameter s), we must impose further conditions on the parameters % and 
B\ here we avoid B — because the curve R(T-l, 0) trivially satisfies periodic 
boundary conditions. In particular, for B =/= 0, if % and B satisfy the 
relationship 



k Vl + B 2 + 2B sin t 

then the curve R(H,B) is 27r/fc periodic in the x variable, for k E N. In 
the sequel, we will use the notation R(B,k) to denote the 2ir/k periodic 
undulary curve with free parameter — 1 < B < 1 and parameter H = H(B) 



fixed according to (3.2) 



All of the figures contained herein were generated with the program GNU Octave, version 
3.4.3, copyright 2011 John W. Eaton, and GNUPLOT, version 4.4 patchlevel 3, copyright 2010 
Thomas Williams, Colin Kelley. 
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Figure 1 . Profile curves for a family of spheres and a catenoid, respectively. 
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Figure 2. ir periodic nodary curves with B = 1.03 and B — 1.1, respectively. 




Figure 3. Families of 2tt periodic undulary curves with selected 
parameter values from B = — .99 to B = 0.99, as indicated. 



4. Stability Of Cylinders With Large Radius 

As seen above, the constant function r(x) = r+, for > 0, is an equilibrium 
of (2.1). Moreover, the constant function r(x) = r+ is associated to the cylinder 



r(rv) with radius r*, which is a stationary solution of the original surface diffusion 
problem ( 1.1 1. In this section, we establish tools for and carry out the investigation 



of nonlinear stability for these equilibria. 

4.1. Preliminary Analysis and Definitions. Throughout this analysis, we con- 
sider an arbitrary > and a 6 \ Z, unless otherwise stated. Focusing on 
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the properties of solutions near r*, we shift our equations, including the shifted 
operator 



G*(p) :=G(p + n) 



1 



d x 



V 1 + pt 



which maps p G E\ D U* to E i where we consider p = r — r*, and is in the regularity 
class by Lemma 2.1 here we take {/* := V — r* := {p — : p G V^}. Now we 
consider the surface diffusion problem shifted by r+, 

(4.1) 



where po 



l/>t(* > a:) = G*(p(* > a)) J i€l,t>0, 
\p(0,a;) = po{x), x € T, 

ro — r*. We say that 

p = p(-, p ) g i, B ) n c( j, Ex) n c(j, f p n £/*) 



is a solution to (4.1 1, with initial data po G F M D [/*, on the interval J C K+ if p 
satishcs (4.1) pointwisc, for t > 0, and p(0) = po- We will investigate the mapping 
properties of G* around in order to gain information about the stability of in 
(pl 



Define the functional 

F*(p)=F i <(p;r ir ) 



(p(x) + r±) 2 dx. 



which corresponds to the volume enclosed by the surface T(p + r+). Considering 
the regularity of F*, it follows from the analyticity of multiplication and integration 
on little-Holder spaces that F* is of class C u from /^(T) to R, cr G M+ \ Z. The 
Frechet derivative of Ft is 



(4.2) F>F*(p) 



(p(x) + r it )h(x)dx 



g £(/i CT (T), 



p e fc, CT (T,R). 



Moreover, it holds that F*(p) is conserved along solutions to (4.1). Indeed, if 
p = p(-, po) is a solution to (4.1 ), then 



d t F*(p(t)) = 2 / (p(t) + r*) Pt {t)dx = 2 d x 



(p(t) + r») 



d x U(p(t)+ rir ) 



dx = 0, 



for t G J{po) \ {0}, where the last equality holds by periodicity. Thus, conservation 
of Ft along the solution p follows by continuity of Ft and convergence of p to the 
initial data po in En. From these properties, it follows that 



(4.3) 



Ml 



{pG^(T):F4p)=F4r;)}, ry 



cr G 



\Z 



is a family of invariant level sets for (4.1 ). The following techniques are motivated 
by results of Prokert [38] and Vondenhoff [45] , whereby one can take advantage of 
invariant manifolds in order to derive stability results. 
First, we introduce the mapping 

F p := p - 



2tt 



p(x)dx, 



which defines a projection on h a (T). We denote by /iq(T) the image F (/i CT (T)), 
which exactly coincides with the zero-mean functions on T in the regularity class 
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h a (T), and we have the topological decomposition 

h a (T) = h5(T) 8 (1 - P )(h lT (T)) SAJ(T)eH. 

In what follows, we equate the constant function [rj(x) = rj\ G (1 — Pq) (/i ct (T)) with 
the value 77 G K, and we denote each simply as 77. 
Consider the operator 

1>(P,P,V) ■= (PqP-P, F*(p)-F*(r?)), 

which maps /V^T) x /i£(T) x K to h£(T) x K and is of class C u , by regularity of 
the mappings F* and P . Notice that $(0, 0, 0) = (0, 0) and, using ( |4.2[ ), 

(4.4) D x $(0,0,0) = (p , 4tit*(1 - ft)) e C isom {h°{T),hZ{T) x K), 

i.e. the Frechet derivative of $ with respect to the first variable, at the origin, is a 
linear isomorphism. Hence, it follows from the implicit function theorem that there 
exist neighborhoods (0, 0) G U = U x U x C h%(T) x R and G U 2 C h a (T) and a 
C u function ip : U —> U2 such that, for all (p, p, rj) G x [/, 

3>(p, p, 77) = (0,0) if and only if p = ^(p, 77). 

Remarks 4.1. We can immediately state the following properties of ip, which follow 
directly from its definition and elucidate the relationship between Pq and ip. 

a) P ip(p,i]) = p for all (p,?/) G U. 

b) Given p G ^>(f7) n it follows that ip(P p, ??) = p . 

c) "0(0,77) = 77, for 77 G t/i. This and the preceding remark follow from the fact 
that F*(t7) is injective when restricted to 77 G (— r*, 00) C E. 

d) It follows from the identity $>(ip(p, 77), p, 77) = (0,0) and differentiating with 
respect to p that £>i$(^(0, rj), 0, r?)Di^(0, 77) - (/i,0) = (0,0). From this 
observation, and the fact that Pi$(?7, 0, 77) = (P ,47r(7\t + 77) ( 1 — Pq)), it 
follows that 

D 1 ip(0,ri)h = h, heho(T), r] eUx. 

e) ip(Uo,r]) C for 77 G E/j.. Hence, ip(-,rj) can be taken as a (local) 
parametrization of M.^ . Moreover, from the preceding remark and the bi- 
jectivity of ip{-, 77) from Uq to M° (1 U2, we can see that M° fl E/2 is a Banach 
manifold over /ip (T) anchored at the point 77 G M . 

f) For (p, 77) G C/, we have the representation 

i>{p, 77) = (P + (1 - P )) ip(p, v)=P+^J 1>(P, ri){x)dx, 

and so we can see that A4° n E/2 can be realized (locally) as the graph of a 
K- valued analytic function over the zero-mean functions p G /io(T). 

g) Although ip(-,rj) depends upon the parameter a, a priori, it follows easily 
from the preceding representation that 

ip(-,r]) : U nh%(T) h a (T), a G M+ \ Z, 

so that ^ preserves the spacial regularity of functions regardless of the regu- 
larity parameter a with which ip was constructed. However, notice that the 
neighborhood Uq will remain intrinsically linked with the parameter which 
was used to construct ip. 
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With the established invariance and local structure of the sets A4°, it follows 



that the dynamics governing solutions to (1.2) reside in the tangent space to the 
manifold M° D 1/2- Hence, if we reduce (1.2) to a local system on M.ZC\ U2, then 



we will have captured all of the dynamics of the problem. Remarks |4.1| ^d) is the 
first observation toward this reduced formulation. In fact, one can make use of the 
properties established in Remarks |4.1| to prove the following, even more general, 
result regarding the properties of the the tangent vectors to A4°. Although we use 



other tools to connect the reduced problem (4.5 1 below with the full problem ( 1.2 ), 



this remark provides good intuition into the nature of these manifolds. 

Remark 4.2. Given (p, 77) G U it follows that Diip{p 1 rf) o P — idj-^,- ^M"i where 
TpAi^ denotes the tangent space to the manifold at the point p. 

4.2. The Reduced Problem. Fix a G (0, 1) and we denote the spaces 

F :=h%(T), fi := h$+«(T), and F„ := (i^OjU' Me (0,1), 

so that F^ = PqE^ for p € [0, 1]. Define the operator 

9*{p,v) = G*(p,V, r *) ■= PoG(i'(p,r]) +7*), 

which is defined for all (p, 77) G U C F x E with p G U D Fi. 

Now we consider the reduced problem for the zero-mean functions 

(p t (t,x) = g+{p(t,x),r 1 ), t>0, 
lp(0,x) =p (x), xeT, 



(4.5) 



where po := -Po r o = Po( r o ~ r *)- One will note that we should insist on ipiPi v)( x ) > 
— r* for all x G T in order to guarantee that G{ip(p, rj)+r^) is well-defined. However, 
we can assume, without loss of generality, that the neighborhood U is chosen small 
enough to ensure this property holds for all (p, 77) G U. 

Remarks 4.3. Throughout most of the analysis that follows, we will treat the pa- 
rameter r\ as a free parameter, although it has a very specific interpretation in 



relation to (2.1). If one is given initial data rp close to r*, then the parameter 77 is 



chosen so that 

F+in) = F*(r ) • 

a) Essentially, this parameter allows for the possibility that the volume enclosed 
by the surface T(ro) differs from that of the cylinder r(r+), thereby allowing 
us to handle non- volume-preserving perturbations r$ of the cylinder r+ . 

b) From a more general viewpoint, one can see that the family {A4° n i/j(U) : 
rj G Ui} forms a dimension 1 foliation of a neighborhood of the positive real 
axis R + C /i CT (T) and the parameter 77 separates the leaves of the foliation. 

For p G (0, 1] and closed intervals J C R + with G J, define the spaces 
E (J) :=BWi-p(J,Eq), 
Ei(J) := BUC'l_^J,E ) n BUCx-^Ex), 

and 

F (J) :=BaCi_„(J,Jb), 

Fi(J) := BUCl_JJ,F ) n BUC 1 - tl (J,F 1 ), 
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within which wc will discuss solutions to the shifted problem (4.1 ) and the reduced 
problem (4.5), respectively. 

In order to connect these two problems, we will make use of the lifting map ip, 
defined in the previous section. To ensure that ip is well-defined on Fi(J), we must 
restrict our attention to functions which map into an appropriate neighborhood 
Uq C Fq of 0. In particular, we assume that Uq is given so that 

tp(-,ri) : Uq C F -> E , V e U u 

is in the regularity class C u and, without loss of generality, we assume that Uq 
is given sufficiently small so that ip and the derivative D\tp are bounded on U = 
UqxUi. More precisely, Uq is chosen sufficiently small so that there exists a constant 
N > for which the inequalities 

(4-6) U(P,r,)\\E <N and \\D^(p, v)\\c(f ,e ) < N 

hold for all (p, r]) G U = Uq x Ui. 

Lemma 4.4. Fix 77 € Ui and J := [0, T] for T > 0. Then 

i>{-,n) :¥ 1 {J)^C{J,Uq) — >Ei(J), with iP(p,r))(t):=ij(p(t),r)). 



Moreover, if po £ F„ and p = p(-,Po) € ^i(J) H C(J,Uq) is a solution to (4.5), 
for some p S [1/2, 1], then p := ip{p,rf) is the unique solution on the interval J to 
(4.1), with initial data po :— i/>(/5o,r/) G E^. 

Proof. First notice that the embeddings 

(4.7) V^J)^ BUC{J,F^)^ BUC{J,F ), p€ [1/2,1], 

follow from [3J Theorem III. 2. 3. 3] and the continuous embedding of little-Holder 
spaces, respectively. 

To see that the mapping property for ip(-,ri) holds, let p £ Fi(J) D C(J,Uo). 
Uniform continuity and differentiability of the function ip(p(-),r]) follows from the 
regularity of ip and p, and compactness of the interval J. Hence we focus on 
demonstrating that ijj(p(-),ri) satisfies the b ound edness conditions for Ei(J). In 
the case p e [1/2, 1), it follows from Remarks 4.1 [f) and (4.6) that, for t E J, 

(4. 8 ) < WpWmj) + * 1- "IIV»(pC0, Eileen 

< IIp|| FiW + T l ~»N, 

and l«at 1 -' , ||^(p(*),»7)||B 1 =0. 

From (4.8) we conclude that ip(p,r]) € BJ7Ci_ M ( J, Ei). Meanwhile, looking at the 
time derivative of ip{p, rf), we note that d t 4>(p(t),i]) = D 1 ip{p(t), rf)dtp{t) and so we 
again make use of ( |4.6| ) to see that 

^-"llWW.jtflk < ||Z> 1 V(p(*),»7)|| J c(j i b, J Bb)t 1 - / *||ftp(*)||*b 
< ^HpIIfxCj) < oo, 
and limt 1 -^!!^^^*),??)!!^ =0. 

Hence, making use of the embedding E\ '-t Eq, we see that ip(p,rj) G Ei(J), as 
desired. Meanwhile, when p = 1 we again get continuity and differentiability from 
the regularity of the mappings p and ip. 
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To see that the second part of the lemma holds, observe by (4.7 1 that po := 
VKpoj ff) € Bp H Hence, by Theorem 2.2 there exists a unique maximal solution 

r(-,po) G C 1 1 _^J(po),E )nC 1 _ f2 (J( Po ),E 1 ) 

to ( |4.1| on some maximal interval of existence J(po) = [0, t + (po)). Now, define 
p(-) := ip(p(-),r}) as indicated and it suffices to show that p t (t) — G+(p(t)) for t £ 
J := (0, T], since this will imply that p(t) = r(t,pa) by uniqueness and maximality 
of the solution r(-,po). Proceeding, let t £ J and consider the auxiliary problem 

f 7 (r) = G*( 7 (r)), forre[0,e], 
17(0)=P(*), 



2.2 



pro- 



which has a unique solution 7 £ C ,1 ([0, e], i?o) H C([0, e], Pi) by Theorem 
vided we choose e > sufficiently small for the particular value p{t) £ E\. Notice, 
by the regularity of 7, we have 

7(0) = G*( 7 (0)) = CM*))- 

Further, note that pit) £ VV^ +Q , from which we conclude that j(t) £ M^ +a and 
by Remarks 4.1 we have the representation 7(7-) = i/j(Po'y(T),rf), r £ [0,e]. Finally, 
we see that 



GMt)) = i(o) = 9 T mp i(T),v)) 



T=0 



D 1 ^(P o7 (0), 77)P o7 (0) 



(4.9) 



= D^iP pit), v )P Q GApit)) = I>iV(p(*),»7)e*(p(t),»?) 

= a t (V(p(*),»j)) = p t (<) J 

which concludes the proof. □ 

We also get the following results which further illuminate the relationship be- 
tween the mappings G* and G*, and explicitly connects the equilibria of the two 
problems (4.1) and (4.5). 

Lemma 4.5. For any p £ A4f. +a fl t/2, it follows that 

(4.10) GAp) = D^iPop, r))PoG*ip), 
and 

(4.11) DG^(p)h = D 2 1 i;iP p,ri)[P h,P G it ip)}+D 1 rP(P p,r 1 )P DG it (p)h, 
for h £ Ei. 

Proof. The first claim was justified in the proof of Lemma[4~4| above and is expressed 



in (4.9). Meanwhile, the second claim follows immediately by differentiation. □ 



Proposition 4.6. If [p, rj) £ U , then (p, 77) is an equilibrium of (4.5) if and only 
if i/j(p,ri) is an equilibrium of (4.1), i.e. 

<=* G*(V(p,»7))=0. 
Moreover, if Q+ip, 77) = 0, then it follows that 

(4.12) DG^ip,i 1 ))h = D 1 yjip, V )P DG4TPip,r ] ))h, h £ E u 
and 

(4.13) DG^ip,r)))D 1 ^ip,r ] )h = D 1 iPip,r))D 1 g*(p,r))h, h £ F x . 
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Proof. The first claim follows from the definition of and (4.10), while (4.12) is a 
consequence of (4.11) and (4.13) follows from (4.11) and the chain rule: 

DG^{p,r)))D^{p,ri)l = DxiP(p,r))P DG^{p,r]))Dx^{p,r])h 

= D x il>{fi,rj)DxQ i ,{p,fi)h. □ 

4.3. Mapping Properties of DxG*{0, rj). Notice that the points (0,7?) S U are 
equilibria of (4.5), and they correspond to the cylinders T(r+ + 77). We are inter- 
ested in the spectral properties of the linearization of Q+ about these equilibria. In 
particular, we compute the Frcchet derivative 

DxG*(0, r{)h = P DG*(V(0, »?)) £>i^(0, ri)h = P DG*{rj) D^{0, rj)h , 

for h € F% . Hence, by Remarks |4.1[ d) we derive the formula 

(4-14) DM0, ri) = P DG±(r))\ Fi = DG^)\ Fi , 

where the last equality is verified by application of the divergence theorem to the 
linearization 



(4.15) 



DGM = -dl 



(r* + v) 2 



Utilizing the Fourier series representation of functions in /i CT (T), c.f. [291 Proposi- 
tions 1.2 and 1.3], we find the eigenvalues of this linearized operator. In particular, 
for he Ex, 



{\-DG*( V ))h= A + 3j 



1 



(4.16) 



MDGM) 



E 

feez 
k 2 



(r* 
A - k 2 
1 



v) 2 

1 



- k- 



(r* + r]) 2 



(r* + ?7) 2 

fc 2 I : fc e Z 



^ ft-(fc)e fe 
h(k)e k 



Noting that the embedding i?i is compact, it follows that the resolvent 

-R(A) := (A — DGi,^))^ 1 is a compact operator, A in the resolvent set p(DG+(r])). It 
follows from classic theory of linear operators that the spectrum a(DG*(rjj) consists 
entirely of isolated eigenvalues of finite multiplicity, see Kato [551 Theorem III. 6. 29] 
for instance. Hence, o- p (DG*(r))) = a{DG i ,{rj)). 

Remark 4.7. If + 77 > 1, then a{DG i< {rj)) C (— 00, 0], however the spectrum will 
always contain 0. The presence of this eigenvalue can be seen as a consequence 
of the fact that the equilibria r* + 77 are not isolated in the space E\. Hence, by 
passing to the operator Q+, which acts on an open subset of the zero-mean functions 
F\, we eliminate the nontrivial equilibria (since the only constant function in Fx 
is the zero function) and thereby eliminate the zero eigenvalue. In particular, one 
easily computes that 



(4.17) a(DxG^0,v)) = \ k 



1 



(r* + rj)'' 



- k' 



k e Z\ {0} 



v e Ux 



Before we return to the problem ( 1.2 ), we state the following maximal regularity 
result for the linearization DxG*(0,r]). For this result, we define the exponentially 
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weighted maximal regularity spaces 

Wj(R+, u>) := {f : (0,oo) ^ F | [t ^ e ut f(t)\ G F i (R+)}, w € R, j = 0, 1, 
which are Banach spaces when equipped with the norms ||w||]j'.(R + ,w) := ||c w 'u||f .(]r + ). 

Theorem 4.8. Suppose r* > 1 and ji G (0,1]. There exist nonzero positive con- 
stants S = 5(r+) and uj — u>(r+, 5) such that 

F (R+,cj), F 1 (R +! a;; 

is a pair o/ maximal regularity for — DiG+(0, rj), given any r\ G (—5,8). I.e. the 
property 

(d t - D&frn),-/) G C isom (F 1 (R+,a;),F (R + ,a;) x h^ +a (T)) , 
Zio/rfs uniformly for rj G (—6, S). 

Proof. Fix <5 > so that (—5, 6) C C7j. n (1 — r*, oo). Following the notation and defi- 



nitions of [29], it is clear from the representation (4.15 1 that —DG+(rj) is a uniformly 
elliptic operator from which we see, by [221 Theorem 4.4], that DG+(rj) generates an 
analytic semigroup on h a (T,C) with domain h 4+a (T,C). Since /iq(T, C) inherits 
the topology of h a (T,C) and the projection Po commutes with DG*(rj), the anal- 
ogous resolvent estimates hold for D\Q±(Q, rj) and so we conclude that DxQ^Q.rj) 



generates an analytic semigroup on hfj(T, C) with domain /iq +q (T,C). Moreover, 



from (4.17) we see that type(D\Q ic (Q,rj)) < for all rj G (—6,6), where type(B) 
denotes the spectral type of the semigroup generator B. In particular, it follows 
that 

type(D 1 g,(0, V )) < <0, rj G (-*,*)• 

Now, choose w G ^0, ^"-ap an< ^ ^he remainder of the result follows from [2, 
Theorem III. 3. 4.1 and Remarks 3.4.2(b)] and the restriction of maximal regularity 
to the subspaces /iq(T). Notice, the characterization 7Fi(w) = /ig /i+a (T) follows 
from stability of little-Holder spaces under continuous interpolation. □ 

4.4. Exponential Stability of Cylinders With Radius r* > 1. Our main re- 



sult regarding stability of cylinders in the axisymmetric surface diffusion flow (1.2) 
establishes exponential asymptotic stability of the family of cylinders under pertur- 
bations which maintain the prescribed periodic boundary conditions and symmetry 
about the same axis of rotation (which we are taking to be the x-axis in our set- 
ting). One feature of our result that we point out is the fact that it establishes a 
form of stability which allows for perturbations which are not volume-preserving. 
In particular, we refer to asymptotic stability of the cylinder r(r*) by which we 



mean that small perturbations of r(r*) will have global solutions to (1.2) which 
converge to a cylinder r(r* + rj), where r* =/= r* + r\ in general. 

Theorem 4.9 (Exponential Stability). Fix a G (0, 1), ji G [1/2, 1], so that Aji + a £ 
Z, and r* > 1. There exist nonzero positive constants e = ^(r*), 6 = 5(r+) and 



u) = u)(r*,8), such that problem (2.1) with initial data r$ £ B/^+a (r*, e) has a 
unique global solution 

r(-,r ) G C 1 1 _^(R +1 h a (T))nC^^(R + ,h i+a (T)), 
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and there exists rj = r](ro) G (— <5, 5) and M = M(a) > for which the bound 
t^Writ, r Q ) - (r* + rf)\\ h A +a + \\r(t,r ) - (r* + v)\\h^+° < e~ ut M\\r - r^+o 
holds uniformly for t > 0. 



Proof, (i) Let 5, w > be the constants given by Theorem 4.8 and consider the 
operator 

£(P> Po, ??) := (dtp - G*{p, ?7), TP - Po) , 

acting on U := (Fi(R+,u;) n C(R+, J7 )) x f^o H F M ) x Ui which is open in the 
Banach space Fi(R + , w) x x K. 

Considering the space into which /C maps, first notice that 

7 :F 1 (R + ,u;)^(F ,F 1 )2 )OO 

follows from [9l Lemma 2.2(a)]. Hence, 7/5 € F M and dt maps Fi(M+,w) into 
Fo(K+,a;) by definition of the spaces BUC\_AJ,E). To see that <?*(•, 77) maps 
U into F (M+,w), choose p G U and notice that p(t) G U n /ig +Q (T), for f > 0, 
from the embeddings (4.7 1. Utilizing the explicit quasilinear representation of the 
operator G, as given by (2.2)-(2.3), whereby 

&(p(t), ??) = Po ( - -4(V'(pW, r?) + r,) (V(p(*), »?) + r*) + f(ip(p(t),v) + r*)) , 

one will easily conclude the desired mapping property for the operator For 
instance, we have seen that A(p)p — bi(p)d^p + ^(pjo^p, where the functions 6, 
only depend on p, p^ and i = 1,2. Hence, it follows that 

e w V-" ^(^(p(t), r?) + r*) (^(p(*),7|) + r*) 

<e-V-^||^V(p(i),r7)|| £o ||&i(V'(pW^) + ^)L 
+ rt 1 -" 11^(^)^)11^ Hfe^C*)^) +r*)|| Eo , 

for t > 0. From here, we take advantage of the boundedness of "i/)(p(t),r)) in the 
topology of Fx/2, m conjunction with the explicit formulas for fy, in order to bound 
the terms \\bi(ip(p(t),i]) + r * ) 1 1 ^ , uniformly in t. Meanwhile, the representation 
given by Remarks |4.1[d) and the fact that p G Fi(K + ,a;) yield the bounds 



e^t 1 -^\^{m,v)\\E =e^t 1 -^mU < l|e""p1k( K+ ), k = 1,...,4. 

Analogous methods work for the remaining terms of the function G*(ip(p(t), T])), 
since we can always isolate an element of the form d^t/j(p(t),r]), and bound the 
remaining elements using boundedness in Fu2- We conclude the result by noting 
that the linear projection Po adds no complexity to acquiring the necessary bounds. 

With the establishment of the spaces into which the operator JC maps, we move 
on with our analysis. Regarding the regularity of K,, it can be shown that G+ is C u 
via substitution operators and the derivative dt and the trace operator 7 are linear. 
Hence, it follows that 

/CgC"(u,F (K + ,w)x V). 
Meanwhile, notice that /C(0,0,0) = (0,0) and 

Di/C(0,0,0)= (ft-Di&M.T) GAaom(Fi(R + ,w),F (R+,a;)xF M ), 
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by Theorem |4.8| Hence, we conclude from the implicit function theorem that there 
exists an open neighborhood G U C-F^xl and a C u mapping n : U — > Fi(R+,w) 
such that 

£( K (/5o, V), Po, V) = (0, 0) for all (p , »y) € l>- 

In particular, n(po, tj) is a global solution to ( |4.5[ ) with parameter 77 and initial data 
Po G Pa, where we assume, without loss of generality, that U CU. 

(ii) Choose e > so that for every r £ Bp (r*, e), there exists 77 G (— r - *, 00) for 
which 

(Por ,r)) G E7 and F*(ro - r^jr*) = F*(?7;r*). 

The existence of such a constant e is guaranteed by the continuity of Po and P*, 
injectivity of F*(t7; r*) for 77 G (—7**, 00) and the fact that Po r * = 0. 

Let ro G 1^(^,6) and fix 77 = rjifo) as mentioned so that F+(ro — rv;r*) = 
.F*(?7;r*). Define the function 

(4.18) r :=-0(K(P o ro,??),?7) + r-*, 

where tp(K(Poro,rf),rj)(t) := ^(k(Po?-o, i7)(i), 77), and we will demonstrate that r 
satisfies the desired properties claimed in the theorem. 



To see that r is the unique global solution to (2.1 ) with initial data ro, first fix 
T > and consider the interval J := [0, T]. By the choice of e > we know that 
(Poro,?y) € U and so it follows from part (i) above that K(P r ,ri) G Fi(R + ,w). 



From this we see that k^PoTq,?]) G Fi(J) is a solution to (4.5) with initial data 
Po^o G F m . Thus it follows, by Lemma 4.4 that r G Ei(J) is the solution on J to 



the problem (2.1) with initial data 

^(P r , 77) + r* = ?K-Po(n) - f+),77) + r* = r , 



where we use Remarks 



4.1 



'b) and the fact that ro — r* G -M^ 1-1 "". The claim now 
follows by the fact that J' > was arbitrary and by definition of the Frechet spaces 

Now, to see that r satisfies the exponential bounds in the second part of the claim, 



first notice that k(0, 77) = for 77 G U\. Then, by definition of r, Remarks 4.1 and 
application of the mean value theorem, we see that the expression 

r(t)-(r* + 77) = ip(n(P r , v)(t),v) - V = ^(^(Po r o, »?)(*), rj) - ^(k(0, 77) (t), 77) 

= (fb + (1 - fo)) (V(K(iVo, »7)(*),»7) - lK«(0, »»)(<), »?)) 

= K(P r ,rj)(t) + ~ J (ii(K(P r ,r))(t,x),r)) - ip(n{0,r})(t,x),r})^ 



dx 



K{Por Q ,ri){t) + — / / D 1 tp(TK(P Q r ,r])(t),T])K(P ro,f])(t,x)dTdx, 



T JO 



2tt 

holds for all t > 0. Notice that 

e"V-^|| K (Poro,77)(<)|| Fl < ||/s(i%r , »?)|| Fl (B 

and 

sup 1 1 e^(P r , 77) (<) 1 1 Fm 
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are finite quantities by the fact that /c(Po r Oj?7) G Fi(R+,u;) and the embedding 
(4.7). We note that the reference for ( |4.7| ) does not explicitly include the un- 
bounded interval J = R+, however the methods of the proof extend to this un- 
bounded case with little trouble. Meanwhile, the remaining term in r(t) — (r* + rf) 
above is scalar- valued, so we bound Diip(TK,(Por , T])(t), rj)K,(P Q r 0l f])(t) in the C(T)- 
topology, which are then bounded in the /i cr (T)-topology, for any a G M. + \ Z. In 
particular, observe that, by (4.6), 

sup ||X>iV(p,»»)K(Poro,»7)(t)IU- < N\\K(P r ,r])(t)\\ hs , t > 0, 
and we conclude that the bounds 

(4.19) e"H 1 -»\\r(t)-(r* + V )\\ El < (l + Cl iv) \\K(P r , 
and 

(4.20) e^\\r{t)-{r,+ r,)\\ E ^ < (c 2 + c 3 iv) \\ K (P r , »?)||f 1 (r + , (U ), 

hold uniformly for t > 0. Here the constant c\ comes from the embedding F\ Fq, 
and the constants C2 and C3 come from the embeddings (4.7). Finally, by the 
regularity of k, we may assume that U was chosen sufficiently small to ensure that 
Dik is uniformly bounded from U into Fi(M + ,w). Recalling that k(Q,t)) = 0, it 
follows that 



(4.21) 



\K{P r ,r])\\v 1 (M + ,u 1 ) 



< 



\D 1 K,(TP r ,r])Poro\\ ¥i 



dr 



< M\\P r \\ Fu < M\\ro-rJ\ Eu 



where M := \\P \\ su P(g,^ e & ll-Pl^P; ??)I U(j- m ,f 
and the inequalities (14.196— (|4.20b. 



The claim now follows from 

□ 



5. Instability of Cylinders with Radius < < 1 



In this section we turn our attention to the stability of cylinders with small 



radius. Again, taking advantage of the reduced problem (4.5) and the connection 



we have established between it and the original problem ( 1.2 1, we proceed with the 
following result regarding instability of cylinders with radius < r+ < 1, in the 
setting of F^. Because of the nature of the instability result, differences in volume 
between the initial data ro and the cylinder r* are not a factor in the following 
argument. In light of this, we will assume that the parameter 77, associated with 



the reduced problem (4.5), is simply taken to be zero for this proof. 



Theorem 5.1. Let r+ E (0, 1) and fj, 6 [1/2, 1] be fixed with + a ^ Z. Then the 
equilibrium of (4.5) is unstable in the topology of /i M+Q (T) for initial values in 
T). 



Proof, (i) Let r± g (0,1) be fixed, and let L :— 0x0^(0,0) be the linearization of 
G* at p — 0. We can restate the evolution equation (4.5 1 in the following equivalent 
form 



(5.1) 



Pt~Lp = g(p), 
P(0) = Po, 



t > 
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where g(p) := G*(p, 0) — Lp. Using the quasilinear structure of [p i— > 0)] it is 
not difficult to see that for every /3 > there exists a number e = £ o(/?) > such 
that 

(5-2) \\g(p)\\ Fo < p\\p\\ Fl , peB Fi ,(0 > e o )nF 1> 

where we will be assuming throughout that p€ [To, to guarantee that G+ip, 0), and 
subsequently g(p), is defined, ft follows from ( |4.f 7 1 that 

o-(i) n [Rez > 0] ^ 0, 

and we may choose numbers u), 7 > such that 

[uj - 7 < Re 2 < w + 7] n a(L) = and 0+ := [Re z > uj + 7] n a(L) ^ , 

i.e. the strip [u> — 7 < Re z < uj + 7] does not intersect a(L) and there is at least 
one point of <r(L) to the right of the line [Rez = w + 7]. 

We define P+ to be the spectral projection, in Fq, with respect to the spectral 
set er + , and let P_ := 1 — P + . Then P + (Fo) is finite dimensional and the topological 
decomposition 

F = P + (F ) ® P_(F ) 

reduces L, so that L = L + fflL_, where L± is the part of L in P±{F ), respec- 
tively, with the domains D(L±) — P±(F\). Moreover, P± decomposes F\ by the 
embedding F\ Fq, and, without loss of generality, we can take the norm on F\ 
so that 

We note that 

cr(L_) C [Rcz < uj - 7], cr(L + ) = (T + C [Rez > uj + 7]. 
This implies that there is a constant Mo > 1 such that 
lle^'P-IU^Moe^)', 
IK L+t F + || £(fo) <M e-(«+^, t>0, 

where {e L_t : t > 0} is the analytic semigroup in P_(Fo) generated by L_ and 
{e L+t : t e R} is the group in P + (F ) generated by the bounded operator L + . 

From ( |4l4| - pT5l ) and [Ml Theorem 5.2] one sees that (F (J),Fi(J)) is a pair 
of maximal regularity for — L and it is easy to see that — L_ inherits the property 
of maximal regularity. In particular, the pair (P_(F (J)), P_(Fi(J))) is a pair of 
maximal regularity for — L_. In fact, since type(—uj + L_) < —7 < we see that 
(P_(F (R + )),P-(Fi(R + ))) is a pair of maximal regularity for (uj - LJ). This, in 
turn, implies the a priori estimate 

(5.4) \\e-^w\\ ¥l[jT) < Mx(|Klk + ||e- w */l|F (^)) 

for Jt := [0, T], any T £ (0, 00) (or Jt = M+ for T = 00), with a universal constant 
M\ > 0, where w is a solution of the linear Cauchy problem 

I w — L_w = /, 
\w(0) = w , 



with (J,wo) E (C((0,T),P-F ),P-U 
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(ii) By way of contradiction, suppose that the equilibrium is stable for (4.5). 
Then for every e > there exists a number d > such that (5.1) admits for each 
po G Bj^(0, S) a global solution 

p = p(;Po) G C 1 1 _ tl (R + ,F )nC 1 ^ l (R + ,F 1 )nC(R+,Uo), 
which satisfies 

(5.5) \\m\\F. < e, t>0. 

We can assume without loss of generality that /? and e are chosen such that 

(5.6) 2Co(M + M l7 )/3 < 7 and e < e (/3), 

where Co := max{||P_||£( Fo ), ||P+||,c(_f )}- As P + (F ) is finite dimensional, we may 
also assume that 

\\P + v\\ F „ = \\P+v\\ Fo , vGF , u€{fi,l}, 

where we also use the fact that P+Po C D(L n ) for every n E N, c.f. [301 Proposition 
A.1.2]. 

CLAIM 1: For any initial value po G M F (0,<5), P + p admits the representation 



(5.7) 



e L + {t ~ s) P+g{p~{s)) ds t>0. 



For this we first establish that, for po e M F (0, S), 

e-^eSCi-^R+^i) := juGC((0,oo),P0 

First notice that the mapping property 

g:¥ 1 (J T )nC(J T ,U )^¥ (J T ), 

which follows in the same way as the mapping property derived for in the proof 
of Theorem 4.9 above, together with the inequalities (5.2) and (5.4) yield 

(5.8) 



sup t 1 ^u{t)\\ Fl < oo 



< T < oo, 



P-P\\b 1 _„(j t ,F!) 



< 



M 1 (||P_po||F M +Co^||e--*P + p|| Bl _ M(jT>Fl) +Co^||e- wi P_p|| Bl _ M(Jr>Fl) ) 



for any < T < oo. Due to (5.6), we have MiCqP < 1/2 and can further conclude 

(5.9) We-^P-pU^j^) < 2A/ 1 (||P_p || F(l +Co/3|| e - wt P + p|| Bl _ M(jT , Fl) ). 

It follows from ( |5.5[ ) that 

t l -^\\e- Mt P + p{t)\\ Fl < ^"^e— f C ||p(*)]|^ < C C ie 



where Ci := supji 1 M e : t > 0} < oo. Inserting this result into (5.9) yields 



(5.10) 



\\e- wt p\\B^(j T , Fl ) < 2M 1 ||P_po||f (J + (2M 1 C 0/ 3 + l)C C x e < C 2 



for any < T < oo. However, since T is arbitrary and (5.10) is independent of 
T we conclude that e'^p G PCi_ /i (M+, Pi), for any initial value p G 18^(0,(5). 
Next we note that, for s > t, by (5.3) 



(5.11) 



\\e L ^P + g(p(s))\\ Fo < M C f3e^^\\p(s)\\ 



<M Coi8e ut e*i , ( t -)^- 1 ||e- w 'p||B 1 . 
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which shows that the integral in (5.7 1 exists for any t > 0, with convergence in F±. 
Moreover, 



(5.12) 



e L + (ts) p+g( ~ {s))ds 



< 



e^M C C 3 p\\e-^p\\ Bl _ 



F 



where C3 := sup { J t °° e 7 '* s 's AI 1 ds : t > 0} < 00. Noting that w = P +j 5 solves 
the Cauchy problem 

\w- L + w = P+g(p), 
\w(0) = P+p , 

it follows from the variation of parameters formula that, for t > and r > 0, 

P+ffl = e L +^P + p(r) + j e L +^P + g(p(s)) ds. 

Since this representation holds for any r > 0, the claim follows from (5.3) and (5.5) 
by sending r to 00. 

CLAIM 2: If p € B F(1 (0,5) and \\p(t, p )|| Ff< < £ for all t > 0, then it must 
hold that 

WP+PoH <2M M 1 Cs\\P-p \\ F .. 
From ( [5J ]) and (|5.11|) follows 

(5.13) 



e wt -P+/5||s 1 _ M (K + ,F ) 



< 



MqCq/3 
1 



t P+/5|U 1 _,( K+> F I ) + \\e- ut P-p\\ 



-,Fi)) 



where we have used the fact that sup t>0 [t 1 M J.°° e 7 *-* s 's^ 1 ds} < 1/7. Adding 
the estimates in (5.8) and (5.13) and employing (|5.6|) yields 



le-^pllB,. 



, Fl) <2M 1 ||P_p || Ffi 



(5.14) 

The representation (5.7) in conjunction with (5.12) and (5.14) then implies 

(5.15) \\P+Po\\f, < M CoC 3 /3\\e'^p\\ Bl _ AR+ . Fl) < M C 3 \\P.po\\ Flx , 

where the last inequality follows from the fact that 2CqM\P < 1. We have thus 
demonstrated the claim. 

Notice that the preceding claim contradicts the stability assumption. In par- 
ticular, if po £ Mp (0,5) is chosen such that P-p~o = 0, then it must hold that 
P+Po = 0, and hence pa — 0, which contradicts the assumption of stability for 
arbitrary po £ B F (0, 5). □ 

We also state the following corollary, which establishes instability of small cylin- 
ders for the original problem (1.2). The corollary is easily proved by use of the 
projection Pq, which serves as a connection between the problems (1.2) and (4.5), 
and application of the result established in the theorem above. In fact, the same 
techniques used to prove Theorem 5.1 can also be employed to prove the corollary 
directly. 

Corollary 5.2. Let r* £ (0, 1) and p £ [1/2, 1] be fixed with Ap + a ^ Z. Then the 
equilibrium r* of (1.2) is unstable in the topology of h^ l+a (T) for initial values in 
h^ +a (T). 
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6. Bifurcation Results 



In this section we turn our attention back to the general equilibria of (1.2 1. In 
particular, we are interested in the interactions between the family of cylinders and 
the family of unduloids. We have already seen that the radius r* = 1 plays a critical 
role in the dynamics of the cylinders. The change of stability for cylinders above 
and below this critical radius suggests that there is a bifurcation at r* = 1. Indeed, 
we will confirm this bifurcation, using results of Crandall and Rabinowitz [13] , and 
investigate properties of the bifurcation. Herein we take the parameter A := 1/r* 
as our bifurcation parameter, r* > 0. 

With the tools and reductions developed in Section |4j we see that it suffices to 
study the bifurcation equation 

(6.1) G(p,\):=G*{p,0) = FbG(il)(ft)+r i ,) = 0, A = 1/r*, 

in the setting of (p, A) 6 Fi x (0, oo), where we use ip(p) '■= ip(p,0) to economize 
notation. Recalling the explicit characterization ( 4.17| ), we note that the eigenvalues 



of -DiC7*(0,0) all have multiplicity two in the setting of F%, regardless of the value 
of the parameter r*. From this observation we see that the techniques of Crandall 
and Rabinowitz |13j . where the authors derive results for operators with simple 
eigenvalues, are not directly applicable in this setting. We may choose at this point 
to employ more general bifurcation results for high dimensional kernels, such as the 
results contained in Kielhofer [37J Section 1.19], or we can simplify the setting in 
which we are working in order to make accessible the results of |13j . 

Whether we choose to simplify our current setting or use the higher dimensional 
bifurcation results, we can make good use of the following observation. Due to the 



periodicity enforced in the problem, the set of equilibria of (1.2) is invariant under 



shifts along the axis of rotation. More precisely, recalling the translation operators 



T a as discussed in the proof of Theorem 2.4 one can easily verify that G(T a r) = if 



and only if G(r) — 0, a <E K. Obviously, this invariance carries over to the reduced 



problem (4.5) and subsequently to the bifurcation equation (6.1). 

One can take advantage of this shift invariance of equilibria in the context of 
bifurcation with high dimensional kernels by constructing a two dimensional bifur- 
cation parameter A = (1/r*, a) and eventually observe two dimensional bifurcating 
surfaces of equilibria in Fx, c.f. 27, Theorem 1.19.2 and Remarks 1.19.3]. On the 
other hand, we will make use of this invariance to simplify the setting in which 
we are looking for equilibria and make accessible the methods of Crandall and 
Rabinowitz for operators with simple eigenvalues. The specific simplification that 
we apply to our setting has also been employed by Escher and Matioc [TH] and is 
supported by the following proposition which allows us to consider the class 

of functions which are even, i.e. symmetric about [a: = 0], and h^ +a regular. 



Proposition 6.1. For every equilibrium p of (4.5), there exists x$ = xo(p) G T 
for which the translation T XQ p is in the space F\ e := /ig*"(T) of even functions 



on T in the class F%. I.e. up to translations on T, all equilibria of (4.5) are even 
functions. 



Proof. From Remarks |3.2| and Proposition |4.6[ we know that p must correspond 
with the projection of an undulary curve R(^H.,B), modulo translations along the 
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x-axis. Choose xq £ T so that T Xo p = PoR(- ;W,B) and one readily verifies that 
R(- ;H, B) is symmetric about s = tt/2H. The claim now follows from x(ir/2'H) — 
0. □ 

From this observation, we see that there is no loss of generality if we focus our 
bifurcation analysis on the setting of p £ F\ e . One benefit of working in this setting 
is that we have the Fourier series representation 

p(x) — ^^a k cos(kx), {a k } C M for all p£F lyC . 

k>l 

We are now prepared to prove our first bifurcation result. 

Theorem 6.2 (Bifurcation of Reduced Problem). For every I £ N, (0,£) £ 
/iq~^ q (T) x (0, oo) is a bifurcation point for the equation (6.1). In particular, there 
exists a positive constant Si > and a nontrivial analytic curve 

(6.2) {(&(«), X e (s)) £ hl+ a xR:s£ (S t ,S t ), (p e (0),X e (0)) = (0,£)} , 

such that 

Q(pt{s), Xt(s)) = for all sg(-fc^), 



and all solutions of (6.1) in a neighborhood of (0,£) are either a trivial solution 
(0, A) or an element of the nontrivial curve (6.2). Moreover, if X £ (0, oo)\N, then 
(0, A) is not a bifurcation point for (6.1). 

Proof. We first note that bifurcation can only occur at points (0, A) for which 
F>iG{0, A) is not bijective. We can see from ( |4.14[ )- pT5 1 that 



(6.3) 



F>iG{0, X) 



-dl(X* + dl 



Fi t , 



which is realized as a Fourier multiplier with the symbol 

and we see that the operator is bijective whenever A £ (0, oo) \N. Hence, it follows 
that bifurcation can only occur at points of the form (0,^), I £ N. 

Now fix £ £ N and we proceed to verify that (0, £) is indeed a bifurcation point 
for |6T| ). By comp_actness of the resolvent R(X) := (A— £>G*(0)) _1 , A £ p(DG il (0)), 
it follows that DiQ(0,£) is a Fredholm operator of index zero. Further, we see that 

N ( := N{D 1 g{0,£)) = span{cos(&)}, 

Ri := R(DiG(Q, £)) = span {cos(fcr) : k > 1, k ^ £} , 

where N(B) and R{B) denote the kernel and the range, respectively, of the operator 
B. Since h a (T) °-> ^(T), we can borrow the L2-inner product to realize Ni as a 
topological complement to Ri as subspaces of Fi >e . Meanwhile, following from ( 6.3 ), 
we compute the mixed derivative 



(6.4) 

Now take vq 



F lt . 



cos(£ •) £ Ni and observe that 

IhDiQQ, i)v = 2£ 3 cos{£ •) £ R { , 
from which the result follows by [131 Theorem 1.7], or [57J Theorem 1.5.1]. 



□ 
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Remark 6.3. Following from the previous result, we are able to track the behavior 
of the so-called critical eigenvalue pte{X) of the linearization DiQ(0, A) about the 
trivial equilibria (0, A). In particular, we choose fJ-e(X) to be the eigenvalue of 
DiQ(0, A) which passes through with non-vanishing speed at A = £, the existence 
of Hz (A) is guaranteed by the bifurcation observed above, c.f. [271 Section 1.6 and 
1.7]. Moreover, employing eigenvalue perturbation techniques, we can also track the 
associated perturbed eigenvalue pe(s) of the linearization DiQ(pe(s), Xi(s)) about 
the nontrivial equilibria. These eigenvalues will play a crucial role in the following 
instability results for the branches of bifurcating equilibria. 

Theorem 6.4. Each of the bifurcations established in Theorem \6.2\ is a subcritical 
pitchfork type bifurcation. More precisely, for all £ € N, we have 

Xt(0) = and A £ (0) < 0, 

where " ' " denotes the derivative with respect to the parameter s. Moreover, it 
holds that the perturbed eigenvalues fj-e(s) are strictly positive for \s\ > chosen 
sufficiently small. 

Proof. Utilizing the methods of [27l Section 1.6 and 1.7], and the techniques de- 
veloped in the previous sections of the paper, one can explicitly verify that the 
bifurcations observed above are indeed subcritical pitchfork bifurcations. To begin, 
for I £ N fixed, recall the subspaces Re, Ne of Fi >e , and io = cos(£-) as defined in 
the proof of Theorem |6 . 2 1 above . First, one computes the second derivative 

Dlg{0,£)[v o ,v ] = -6^ 5 cos(2f) e Re, 

from which we conclude A^(0) = 0. Meanwhile, utilizing the representation |27l 
(1.6.11) and (1.6.9)] and following a considerable amount of computation, one will 
see that 

Mo) = -~<o. 

6 

From these observations, we conclude that the bifurcation observed at the point 
(0, £) is a subcritical pitchfork bifurcation. 

The result for the perturbed eigenvalues now follows from the eigenvalue pertur- 
bation techniques in [571 Section 1.7], see also Amann [TJ Section 27]. □ 

With these bifurcation results established in the setting of the reduced problem, 
we will now go about deriving results for the original problem (1.2 1. Recalling the 



definition of the operator G* from Section [4. 1| we introduce the notation 

G(p,\):=G(p+l/\) = G*(p), for A = 1/r*. 
We are now interested in finding solutions to the bifurcation equation 
(6.5) G{p, A) = 0, {p, A) e h i+a (T) x (0, oo), 



associated with the full problem (1.2). 

We begin analyzing (6.5) by lifting the bifurcation results already established 



for the reduced problem. We make use of the connections established in Section 4.2 



and we also establish the following connection between the eigenvalues of D\Q and 
DG(-,X) at equilibria. 

Proposition 6.5. Suppose Q(p,X) = and /.i ^ 0. Then 
p is an eigenvalue for D\Q{p, A) ■<=>■ p is an eigenvalue for DiG{ip{p), A). 
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Proof, (i) First, suppose that DiQ(p, X)h — ph for some h € F\ \ {0}, and let 
h := Dip(p)h. Then h € E±\ {0}, by injectivity of Dip(p), and it follows from 



( |4.13| that 

D x G{ip{p),X)h = ph. 

We also observe that this assertion is true in case p = 0. 

(ii) Now suppose that DiG(ip(p), X)h = ph for some h £ E\ \ {0}. We conclude 



from ( 4.12[ ) that h £ T^(p).Mo, so that there exists a unique h £ -Fi \ {0} for which 



/i = Dtp{p)h. Then ( |4.13[ ) shows that 

pD^(p)h = Dij){p)D x G{p,\)h, 
and finally, by injectivity of Dip(p), we conclude that ///i = £/(/>, A)/i, as desired. □ 

We are now prepared to prove the main result regarding bifurcation of the orig- 



inal problem (1.2 1 in the setting of h 4+a (T), and instability of the bifurcating 
unduloids. 

Theorem 6.6 (Bifurcation of Full Problem). Fix I £ N. Then: 
a) the set 

(6.6) + 1/M*) : « G C /i 4+Q (T), 



«s an analytic curve of equilibria for the problem (1.2 1 which bifurcates sub- 
critically from the family of cylinders r+ £ (0, oo), at t/ie cylinder r± = l/£. 
b) </iere exists some eg, > so t/iat /or every s £ {—St, 5g) 

Tp{Pe( s )) + !/ A f( s ) = R(B,i), for some B £ (-eg, eg), 
i.e. the family (|6.6[) o/ equilibria are exactly the even presentations of 2ir It- 



periodic undulary curves in some neighborhood of the cylinder = l/i. 
c) the undulary curves R(B,£) are unstable for \B\ > chosen sufficiently 
small. 



Proof, (a) It follows from Proposition 4.6 and Theorem 6.2 that the family 



{(^(pg(s)),\g(s)) : s £ [-6 t ,6t)} C Si x (0,oo) 



consists of solutions to the bifurcation equation (6.5). The regularity of the curve 
follows from the regularity of the bifurcating branch in Theorem |6 . 2 1 and regularity 
of the mapping ip. By definition of the bifurcation function G(-, A), it follows that 



the family (6.6) are indeed equilibria of the original equation ( KL2| which intersect 
the family of cylinders at r* = l/£, when s = 0. Meanwhile, the bifurcation 
parameter A remains unchanged in lifting from the reduced problem to the full 
problem, hence we see that 

X e (0) = and A^(O) < 0, 



from Theorem 6.4 and so we conclude that the given curve bifurcates subcritically. 

(b) By Remarks |4.1| ^f) it follows that tp preserves the symmetry of even functions 
on T, and since pg(s) £ F\, e , it follows that the functions in the family ( |6.6[ ) are 
even on T. Meanwhile, by the characterization of equilibria established in Section[3j 
and the fact that 

i)(p e (o)) + i/\ f .(o) = i/i = R(o,e), 
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it follows that the family (6.6| must coincide with the family of 27r/^-periodic 



undulary curves R{B 1 t), for some continuum of values B € (— ei,si). 



(c) To prove that the unduloids (6.6) are unstable, we mimic the proof of The- 
orem [5TT] in the current setting. In particular, define 

G e {p, s) := G{p + tp(pi(s)), Xt(s)), and 
L e (s) := AG £ (0, s) = DiG^ipiis)), X e (s)), 



acting on functions p € E\. It follows by Theorem |6.4| and Proposition 6.5 that 

<j{L t {s)) n [Rez > 0] ^ 0, 
provided |s| > is chosen sufficiently small. Meanwhile, the operator Ge(-,s) has 



a similar quasilinear structure as G+ and so the analogue to inequality (5.2) is also 
derived for 

grip, s) := G e {p, s) - L e (s)p. 



Utilizing [571 Proposition 1.7.2] and the explicit characterization (4.16) of the spec- 
tra a(DGi,(ji)), we can control the eigenvalues of the perturbed linearization Li{s) , 
so that, for sufficiently small values of \s\ > 0, we can derive the necessary spectral 
gap condition 

[w-7 < Rez < w + 7] Da(Li(s)) = and o+ := [Rez > u + 7] D a(L e (s)) ^ , 

for some 7, u > 0. The remainder of the proof now follows as in the proof of The- 
orem 5.1 with the observation that —Li(s) satisfies maximal regularity properties, 
which follows by uniform ellipticity of Li(s) and an argument similar to the proof 
of Claim 1 in Lemma EH Q 



Remark 6.7. The instability result Theorem 6.6 c) can also be proved in the reduced 
setting for the functions pe{s), again provided that |s| > is sufficiently small. 
However, it is more difficult to prove that the associated linearization Li(s) satisfies 
maximal regularity properties due to the presence of the zero-mean projection Pq 
in the definition of the operators in the reduced setting. 
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